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Summary 

Swine acute diarrhea syndrome coronavirus (SADS-CoV), a novel coronavirus, was first 
discovered in southern China in January 2017 and caused a large scale outbreak of fatal 
diarrheal disease in piglets. Here, we conducted a retrospective investigation of 236 samples 
from 45 swine farms with a clinical history of diarrheal disease to evaluate the emergence 
and the distribution of SADS-CoV in pigs in China. Our results suggest that SADS-CoV has 
emerged in China at least since August 2016. Meanwhile, we detected a prevalence of 
SADS-CoV (43.53%), porcine deltacoronavirus (8.83%), porcine epidemic diarrhea virus 
(PEDV) (78.25%), rotavirus (21.77%) and transmissible gastroenteritis virus (0%), and we 
also found the co-infection of SADS-CoV and PEDV occurred most frequently with the rate 
of 17.65%. We screened and obtained two new complete genomes, five N and five S genes of 
SADS-CoV. Phylogenetic analysis based on these sequences revealed that all SADS-CoV 
sequences in this study clustered with previously reported SADS-CoV strains to formed a 
well defined branch that grouped with the bat coronavirus HKU2 strains. This study is the 
first retrospective investigation for SADS-CoV and provides the epidemiological information 
of this new virus in China, which highlights the urgency to develop effective measures to 
control SADS-CoV. 
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11 INTRODUCTION 

Swine acute diarrhea syndrome coronavirus (SADS-CoV) is a newly discovered coronavirus 
which is an enveloped, positive and single-stranded sense RNA virus with a genome size of 
approximately 27kb (Gong et al., 2017; Pan et al., 2017; Zhou et al., 2018). SADS-CoV 
belongs to the family Coronaviridae which contains four genera, Alphacoronavirus, 
Betacoronavirus, Gammacoronavirus and Deltacoronavirus (Woo et al., 2012; Woo et al., 
2010). So far, six coronaviruses have been identified from pigs, which include porcine 
epidemic diarrhea virus (PEDV), porcine respiratory coronavirus (PRCV), SADS-CoV and 
transmissible gastroenteritis virus (TGEV) that all belong to the Alphacoronavirus genus, as 
well as one betacoronavirus, porcine hemagglutinating encephalomyelitis virus (PHEV) and 
one deltacoronavirus, porcine deltacoronavirus (PDCoV) (Lin et al., 2016; Woo et al., 2010; 
Wesley at al., 1991). Among these viruses, SADS-CoV is the most newly discovered 
coronavirus, which has been first reported in 2017 in China and is considered to be an 
HKU2-related coronavirus with a bat-origin (Gong et al., 2017; Zhou et al., 2018). In January 
2017, SADS-CoV was detected in a swine farm and subsequently spread rapidly to three 
other farms in Guangdong Province and caused the fatal swine acute diarrhea syndrome 
(SADS) characterized by the clinical signs with severe, acute diarrhea and rapid weight loss 
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of piglets. The symptoms of SADS-CoV are similar to those that caused by other swine 
enteric coronaviruses such as PDCoV and PEDV, but SADS-CoV is more harmful than these 
viruses because it has led to the death of almost 25,000 piglets in a short time and resulted in 
more significant economic losses (Sun et al., 2016; Dong et ah, 2015; Zhou et al., 2018). So, 
it is urgent to investigate the molecular epidemiology and transmission patterns of 
SADS-CoV for establishing effective controls for this new coronavirus. 

In the present study, we performed the retrospective PCR testing on diarrheal samples 
from 45 swine farms in Guangdong Province to evaluate the emergence and the distribution 
of SADS-CoV in pigs in China. The prevalence and co-infection information of SADS-CoV 
from eleven SADS-CoV-positive farms was provided. The sequences of SADS-CoV, 
including two complete genomes, five nucleocapsid protein (N) genes and five spike protein 
(S) genes, were also identified and characterized to investigate the phylogenetic relationships 
of SADS-CoV. 


2 I MATERIALS AND METHODS 

2.11 Sample collection 

A total of 236 clinical diarrhea samples including feces and intestinal contents were collected 
from piglets and sows from 45 swine farms in the cities Qingyuan and Shaoguan of North 
Guangdong Province between August 2016 and May 2017 in accordance with the 
recommendations of National Standards for Laboratory Animals of the People’s Republic of 
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China (GB149258-2010). Samples were preserved at -80 °Cfrom the time of original receipt 


until use. 

2.2 I Nucleic acid extraction and molecular diagnosis 

Samples were homogenized in phosphate-buffered saline (PBS) (20 % w/v), frozen and 
thawed three times, then centrifuged for 10 min at 10,000xg. Viral nucleic acid was extracted 
following the manufacturer’s recommendations of AxyPrep™ Body Fluid Viral DNA/RNA 
Miniprep Kit (Axygen Scientific, Inc). The virus nucleic acid was stored at -80°C until PCR 
was performed. A pair of primers (forward primer 
5'-GGTCCCTGTGACCGAAGTTTTAG-3', reverse primer 5'- 

GCGTTCTGCGATAAAGCTTAAAACTATTA-3') was designed to detected SADS-CoV 

TM 

based on the conserved N gene of this virus. One step RT-PCR using PrimeScript One Step 
RT-PCR Kit Ver.2 with Dye Plus (Takara, Biotechnology, Dalian, China) was carried out to 
amplify the target fragments by the following thermal profile of 50°C for 30min, 94°C for 
3min, 35 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, an extension at 
72°C for 30 s and a final step of 72°C for 5 min. Four other diarrheal pathogens including 
PEDV, PDCoV, rotavirus (RV), and TGEV from SADS-CoV-positive farms were also 
tested by RT-PCR according to the previously described methods (Mai et al., 2017; Liu et al., 
2015; Stevenson et al., 2013). 
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2.3 I Amplification of the N gene, the S gene and the complete genome of SADS-CoV 


Specific primer pairs based on reported SADS-CoV strains (GenBank accession numbers: 
MF094681-MF094684) were designed for S genes, N genes and complete genome 
amplifications respectively (Table SI). PCR assays were performed with the following 
thermal profile: 95 °Cfor 5 min, 35 cycles of 95 °Cfor 30 s, 50 °Cfor 30 s, and 72 °Cfor 1 min 
15 s, followed by a final 10 min extension at 72 °G The products were purified following the 
manufacturer’s instructions of Gel Band Purification Kit (Omega Bio-tek, USA) and then 
cloned into the PMD19-T vector (Takara, Biotechnology, Dalian, China) and transformed E. 
coli DH5a competent cells. The positive clones were screened out and sent to Beijing 
Genomics Institute (Shenzhen, Guangdong, China) for further sequencing. 

2.4 I Sequence alignment and phylogenetic analysis 

The nucleotide sequences were assembled and aligned using the DNASTAR program 
(DNAStar V7.1, Madison, WI, USA). Phylogenetic trees were constructed using the 
neighbour-joining method in MEGA 7.0 software with bootstrap analysis of 1,000 replicates. 
Percentages of replicate trees in which the associated taxa clustered are shown as nearby 
branches (Chenna et al., 2003; Tamura, Nei & Kumar, 2004; Kumar, Stecher & Tamura, 
2016). 
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3 I RESULTS 


3.11 SADS-CoV detection 

Of 236 clinical diarrhea samples from 45 swine farms, 53 samples collected from eleven 
swine farms between August 2016 and May 2017 were tested positive for SADS-CoV and 
the locations of these eleven farms were indicated in Figure 1. The farms SM, GL, DL, DE, 
TP, WT, DCD and SC were in Qingyuan, and the farms LS, ZW and YX were in Shaoguan. 
From these eleven SADS-CoV-positive farms, a total of 170 archived diarrheal samples 
including 55 feces and 115 intestinal contents were collected to further investigate diarrheal 
pathogens including SADS-CoV, PEDV, RV, PDCoV, and TGEV. The infection rates of 
these five pathogens in order were 43.5% (74/170), 78.2% (133/170), 21.8% (37/170), 8.8% 
(15/170) and 0% (0/170). SADS-CoV was tested positively both in the intestinal and fecal 
samples and the detection rates were 49.6% (57/115) and 30.9 % (17/55), respectively, which 
were both the second highest among the five tested pathogens. PEDV had the highest positive 
rate at 84.3% (97/115) for intestinal and 65.5% (36/55) for fecal samples. TGEV was tested 
negatively in all 170 diarrheal samples (Table S2). 

The cases of individual infection or co-infection for SADS-CoV, PDCoV, PEDV, and 
RV were also detected. More than half the samples were only infected by one of the four 
viruses at a rate of 58.82% (100/170), and the individual infection of PEDV had the highest 
rate of 37.06% (63/170). In the 74 SADS-CoV-positive samples, 28 samples were infected 
alone by SADS-CoV, and the rest 46 samples were infected by SADS-CoV combined with 
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one to three other pathogens. The co-infection of SADS-CoV and PEDV occurred most 
frequently with the rate of 17.65% (30/170) (Table 1). 

The earliest cases with SADS-CoV-positive infection were found in the farms LS, TP 
and ZW on August 2016 according to Table S2. The intestinal contents from these eleven 
farms were tested positive for SADS-CoV and the positive rates ranged from 25% to 77.78%. 
While the fecal samples from two farms ZW and YX were tested negative for SADS-CoV, 
but positive for PDCoV, PEDV and RV infections. The positive rates of SADS-CoV in the 
fecal samples of other nine farms varied from 13.33% to 100%. The farm LS had the highest 
positive rate of 77.8% (7/9) for intestinal samples and the farms DCD and SC had 100% 
positive rates for fecal samples. 

3.2 I Sequence identities and phylogenetic analyses 

The two complete genomes (accession number MG605090 and MG605091), five N genes 
(accession number MG605087, MG605088, MG605089, MG775251 and MG775253) and 
five S genes (accession number MG605084, MG605085, MG605086, MG775250 and 
MG775252) were sequenced and available in GenBank. The lengths of the complete genome, 
N gene and S gene were 27163bp, 1128bp and 3390bp respectively. The two new complete 
sequences, five N genes and five S genes sequences each shared high sequence identities of 
99.8% - 100%, 99.8% - 100% and 99.6% - 100% with previous reported complete genomes, 
N genes and S genes of SADS-CoV. The two full-length sequences of SADS-CoV showed 
95% - 95.1% nucleotide identities with four bat coronavirus HKU2 strains, and the nucleotide 
identities between the N genes of S ADS- CoV studied here and those of bat coronavirus 
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HKU2 strains ranged from 94.2% to 94.3%. The S genes of SADS- CoV in this paper shared 


low nucleotide identities of 79.9% - 80% with the S genes of bat coronavirus HKU2 strains. 
Phylogenetic analyses based on the complete genomes, N genes and S genes indicated that all 
SADS-CoV sequences clustered together to form a well defined branch and group with four 
bat coronavirus HKU2 strains In the phylogenetic trees of the full-length genomes and N 
genes, all SADS-CoV strains were phylogenetically located within the AlphaCoVs, while the 
phylogenetic tree of S genes exhibited that all SADS-CoV sequences, four bat HKU2 and 
two BtRf-AlphaCoV strains separated from other AlpaCoVs sequences and clustered with the 
BeltaCoV sequences (Figures 2, 3, 4). 


4 I DISCUSSION 

In this study, we performed the retrospective investigation of SADS-CoV in 45 pig farms 
from Guangdong Province based on 236 diarrhea samples. Our results showed that the first 
SADS-CoV positive sample was collected in August 2016 from the farm LS with a history of 
diarrhea, as well as from other two farms TP and ZW, which indicates that SADS-CoV has 
emerged in pigs in China at least since August 2016. And this time point is five months 
earlier than the first discovered time reported by our previous study (Zhou et al., 2018). As 
the same time, clinical signs of SADS-CoV during the retrospective investigation included 
sever and acute vomiting and diarrhea, leading to death in piglets that were less than five 
days of age with a mortality rate of around 50%. These clinical presentations were similar to 
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those signs in the large scale outbreak of SADS-CoV reported by Zhou et al. (2018), except 


the mortality rate in piglets later increased to 90%. 

Based on the rates of infection documented in our work, it revealed that PEDV (78.25%) 
was still the primary cause of the porcine diarrhea, which is consistent with previous studies 
that PEDV has been considered to be the major pathogen responsible for the porcine diarrhea 
epidemic in China since 2010 (Sun et al., 2012; Ge et al., 2013; Zhao et al., 20216 ). Except 
PEDV, SADS-CoV had the second high infection rate (43.53%), which exhibits an evidently 
prevailing tendency in pigs. Meanwhile, of the SADS-CoV positive samples, 62.2% were 
co-infections with one to four other viruses, revealing a high prevalence of co-infection in the 
sampled farms. Notably, our results also showed that SADS-CoV and PEDV were 
simultaneously present in all co-infection samples of SADA-CoV. Considering the fact that 
PEDV has caused prior outbreaks of the porcine diarrhea at the pig farm where Zhou et al. 
(2018) later reported the occurrence of SADS-CoV, PEDV seems to be able to contribute to 
the infection of SADS-CoV in pigs. Thus, SADS-CoV may be pathogenic as secondary 
infection following the infection of PEDV, which needs further studies to better understand 
the pathogenesis of this novel coronavirus. 

The phylogenetic relationships of SADS-CoV sequences were also identified in this 
study. The results showed that all SADS-CoV sequences clustered together to form an 
independent branch and separated from other viral sequences in the genus Alphacoronavirus. 
Our results also indicated that both the complete genomes, N genes and S genes of all 
SADS-CoV strains shared the highest nucleotides identifies with those corresponding 
sequences of four bat coronavirus HKU2 strains. In this work, The phylogenetic trees of full 
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length genomes and S genes of SADS-CoV sequences showed that the SADS-CoV branch 


clustered with these four HKU2 strains, which is same to previous results (Gong et al., 2017; 
Pan et al., 2017; Zhou et al., 2018). Besides the genomes and S genes, the tree of N genes in 
our study revealed the identical result too. So far, a total of eight full-length genomes of 
SADS-CoV have been reported in Guangdong Province of China (Gong et al., 2017; Pan et 
al., 2017; Zhou et al., 2018; this study). The two new genomes of SADS-CoV sequences in 
this work shared 100% nucleotides identities with the sequence MF167434 published by 
Gong et al. (2017) and our four previously reported sequences (Zhou et al., 2018), and shared 
99.8% nucleotides identities with the sequence MF370205 studied by Pan et al. (2017). The 
results suggest that these eight SADS-CoV sequences may come from the same origin. Only 
the phylogenetic tree of S genes in our work showed that sequences of the AlphaCoV were 
divided into two sublineages, AlphaCoVl which contained all SADS-CoV sequences and 
AlphaCoV2, clustering together with sequences of the BeltaCoV and the DelatCov, 
respectively. And this result was consistent with the study of Pan et al. (2017). As a newly 
discovered coronavirus, the availability of SADS-CoV sequences data is limited which 
prevents better understandings of the molecular epidemiology of this virus. Meanwhile, being 
a RNA virus, SADS-CoV may mutate rapidly and exhibit high genetic differences 
(Drummond et al., 2003; Kiihnert et al., 2011). So, new more available sequences data are 
warranted to give a deep insight of viral genetic origin, evolution and transmission patterns of 
SADS-CoV. 
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In summary, our retrospective study suggests that SADS-CoV has emerged in pigs in 


China at least since August 2016. The severe clinical symptoms and the strong transmission 
in short term highlights the urgency to develop effective measures to control this new 
discovered virus. Further studies, such as the epidemiology, virology and pathobiology 
should be performed for better understanding of SADS-CoV in China. 
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Figure Legends 

Figure 1. Locations of positive farms for study of SADS-CoV in Guangdong Province, 
China. 

Figure 2. Phylogenetic analysis of the complete genomes of SADS-CoV and reference 
coronavirus species. The tree was constructed using MEGA 7.0 software with 
neighbor-joining methods and 1,000 replicate sets on bootstrap analysis. Two new complete 
genomes studied in this work were indicated with " black solid circles ". 

Figure 3. Phylogenetic analysis of the N genes of SADS-CoV and reference coronavirus 
species. The tree was constructed as per Figure 2 above. Five new sequences of N genes 
studied in this work were indicated with "black solid circles". 
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Figure 4. Phylogenetic analysis of the S genes of SADS-CoV and reference coronavirus 


species. The tree was constructed as per Figure 2 above. Five new sequences of S genes 
studied in this work were indicated with "black solid circles". 

Table 1. Screening information of SADS-CoV, and co-infections of different pathogens 
detected in north Guangdong Province of China between August 2016 and May 2017. 

Table SI. Primers for amplifications for the competed genome, S and N genes of 
SADS-CoV. 

Table S2. Detection of SADS-CoV, PEDV, RV, PDCoV and TGEV in diarrheal samples 
from eleven farms of North Guangdong Province. 
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Table 1 


Month 

s of 

One pathogen positive samples no.(%, 

n=170) 


Two pathogens positive samples 

no.(%, n=170) 

Three pathogens positive samples no.(%, n=170) 

Four pathogens 

positive samples 

no.(%, n=170) 

SADS-CoV 

positive 

farms 

year 

SADS-C 

oV 

PEDV 

RV 

PDCo 

V 

TGE 

V 

SADS-COV+ PEDV+R 

PEDV V 

PEDV+PD 

CoV 

SADS-CoV+PE 

DV+RV 

SADS-CoV+PEDV 

+PDCoV 

PEDV+RV+P 

DCoV 

SADS-CoV+PEDV+R 

V+PDCoV 

(Positive/T 

otal) 

Aug, 

2016 

3 

4 

0 

0 

0 

2 

3 

1 

0 

0 

1 

0 


Sep, 

2016 

0 

6 

0 

2 

0 

0 

1 

1 

0 

0 

0 

0 


Oct, 

2016 

0 

2 

1 

0 

0 

0 

3 

0 

0 

0 

0 

0 


Nov, 

2016 

1 

8 

1 

1 

0 

1 

4 

1 

1 

1 

0 

2 

11/45 

Dee, 

2016 

2 

8 

0 

1 

0 

3 

0 

1 

3 

0 

0 

0 


Jan, 

2017 

2 

3 

2 

0 

0 

10 

3 

0 

2 

0 

0 

0 


Feb, 

2017 

6 

7 

0 

0 

0 

4 

3 

1 

2 

0 

0 

0 



This article is protected by copyright. All rights reserved. 






Mar, 

3 9 0 1 0 4 0 

2017 


Apr, 

10 10 0 0 0 6 0 

2017 


May,2 

1 6 0 0 0 0 1 

017 


Total 


28(16.4 63(37.0 4(2.3 5(2.9 0(0 
7%) 6%) 5%) 4%) %) 


30(17.65%) 


18(10.5 

9%) 


SADS-CoV positive samples : 74(43.53%) 
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o 


4 


0 


0 


0 


0 0 


0 0 


5(2.94%) 12(7.06%) 2(1.18%) 1(0.59%) 2(1.18%) 





HuNan Province 


JiangXi province 


Shaoguan 


Mcizhou 


Qingyuan 

6-WT, 


Zhaoqing 


JTeyang 


(iuang/ltou 


Huizhou 


Shanwci 




Yunfu 


.Hangmen 


Yangjiang 1 


Maoming 


/.hi njiang 


FuJian Province 
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*E 


• SADS-CoV/CN/GDGL/2017 Complete Genome | MG605090 

• SADS-CoV/CN/GDLS/2017 Complete Genome | MG605091 
Porcine enteric alphacoronavirus GDS04 | MF167434.1 
SADS-CoV FarmC | MF094683.1 
SADS-CoV_FarmB | MF094682.1 
SADS-CoV_FarmD | MF094684.1 
SADS-CoVFarmA | MF094681.1 

Rhinolophus bat coronavirus HKU2 CH/GD-01/2017/P2 | MF370205.1 
Bat coronavirus HKU2 | NC009988.1 
JtI Bat coronavirus HKU2 strain HKU2/HK/46/2006 | EF203065.1 
Bat coronavirus HKU2 strain HKU2/HK/298/2006 | EF203066.1 
Bat coronavirus HKU2 strain HKU2/HK/33/2006 | EF203067.1 
■ BtRf-AlphaCoV/YN2012 | KJ473808.1 
86~l BtRf-AlphaCoV/YN2012 | NC028824.1 
Human coronavirus 229E | NC002645.1 

Camel alphacoronavirus isolate Riyadh/Ry141/2015 | NC028752.1 
229E-related bat coronavirus strain BtKY229E-1 | KY073747.1 
Human Coronavirus NL631NC005831.2 

— NL63-related bat coronavirus strain BtKYNL63-9a | NC032107.1 

- Bat coronavirus HKU8 strain AFCD77 | EU420139.1 

— Bat coronavirus 1A | NC010437.1 

- Rousettus bat coronavirus HKU10 | NC018871.1 

-Bat coronavirus CDPHE15/USA/2006 | NC022103.1 

-Scotophilus bat coronavirus 512 | N009657.1 

t Porcine epidemic diarrhea virus strain GDS01 | KM089829.1 
Porcine epidemic diarrhea virus strain CV777 | KT323979.1 
Porcine epidemic diarrhea virus| NC 003436.1 
p PRCVISU-1 | DQ811787.1 
I Transmissible gastroenteritis virus | AJ271965.2 
65 ItgEV Purdue P1151 DQ811788.1 

57 j- PHEV CoV USA-15TOSU1582 | KY419113.1 
59 L- Porcine hemagglutinating encephalomyelitis virus | NC007732.1 
I— Human coronavirus OC43 strain SC2481 | KY983583.1 
58l— Bovine coronavirus | NC003045.1 

- Murine hepatitis virus strain JHM complete genome | AC000192.1 


-c 

I— M. 


54 I— Mouse hepatitis virus strain MHV-A59 Cl2 mutant | NC001846.1 

- Coronavirus Neoromicia/PML-PHEI/RSA/2011 | KC869678.4 

- Middle East respiratory syndrome coronavirus | NC019843.3 

- Bat coronavirus HKU4-1 | NC009019.1 
-Bat coronavirus HKU5-1 | NC009020.1 


60 


56 ~|_r 

55L 


Bat coronavirus HKU9-1 |NC009021.1 NC 009021.1 
Bat SARS coronavirus HKU3-1 | DQ022305.2 
SARS coronavirus | NC004718.3 
Bat SARS-like coronavirus RsSHC014 | KC881005.1 


55 L Bat SARS-like coronavirus WIV1 | KF367457.1 

- Avian infectious bronchitis | NC001451 

- European turkey coronavirus 080385d | KR822424.1 
-Bulbul coronavirus HKU11-7961 FJ376620.1 

|- Porcine deltacoronavirus isolate PDCoV/CHJXNI2/20151 KR131621.1 
STL Porcine coronavirus HKU15 strain HKU15-44 | NC016990.2 


AplhaCoV 


BetaCoV 
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Porcine enteric alphacoronavirus GDS04 | MF167434.1 

• SADS-CoV/CN/GDTP/2017 N gene | MG605088 

• SADS-CoV/CN/GDZW/2017 N gene | MG605089 

• SADS-CoV/CN/GDLS/2017 N gene | MG775253 

• SADS-CoV/CN/GDGL/2017 N gene | MG775251 

• SADS-CoV/CN/GDDL72017 N gene | MG605087 
Rhinolophus bat coronavirus HKU2 CH/GD-01/2017/P2 | MF370205.1 
SADS-CoV FarmD | MF094684.1 
SADS-CoV FarmC | MF094683.1 
SADS-CoV_FarmB | MF094682.1 
SADS-CoV FarmA | MF094681.1 

r Bat coronavirus HKU2 | NC009988.1 
Bat coronavirus HKU2 strain HKU2/HK/298/2006 | EF203066.1 
Bat coronavirus HKU2 strain HKU2/HK/33/2006 | EF203067.1 
Bat coronavirus HKU2 strain HKU2/HK/46/2006 | EF203065.1 
BtRf-AlphaCoV/YN2012 | KJ473808.1 


99 

100 


100 I BtRf-AlphaCoV/YN2012 | NC028824.1 

-Bat coronavirus 1A | NC010437.1 

-Bat coronavirus HKU8 strain AFCD77 | EU420139.1 

-Rousettus bat coronavirus FIKU10 | NC01887.1 

-Bat coronavirus CDPHE15/USA/2006 | NC022103.1 

-Scotophilus bat coronavirus 512 | N009657.1 

■ Porcine epidemic diarrhea virus strain GDS01 | KM089829.1 


100 


j“‘ 

I-Hu 


iC, 


r Porcine epidemic diarrhea virus strain CV777 | KT323979.1 
72 L Porcine epidemic diarrhea virus| NC003436.1 
Human Coronavirus NL63 | NC005831.2 

-NL63-related bat coronavirus strain BtKYNL63-9a | NC032107.1 

-229E-related bat coronavirus strain BtKY229E-1 | KY073747.1 

Camel alphacoronavirus isolate Riyadh/Ry 141/2015 | NC028752.1 
Human coronavirus 229E | NC002645.1 
j- TGEV Purdue P115 | DQ811788.1 
ISo" j PRCV ISU-1 | DQ811787.1 
1001 Transmissible gastroenteritis virus | AJ271965.2 
541 - Bovine coronavirus | NC003045.1 

83 l— Porcine hemagglutinating encephalomyelitis virus | NC007732.1 
PHEV CoV USA-15TOSU1582 | KY419113.1 
Human coronavirus OC43 strain SC2481 | KY983583.1 
Mouse hepatitis virus strain MHV-A59 Cl 2 mutant | NC001846.1 
Murine hepatitis virus strain JHM complete genome | AC000192.1 

-Bat coronavirus HKU9-1 | NC009021.1 NC009021.1 

- Coronavirus Neoromicia/PML-PHEI/RSA/2011 | KC869678.4 


38 

100 


- Middle East respiratory syndrome coronavirus | NC019843.3 
Bat coronavirus HKU5-1 | NC009020.1 
Bat coronavirus HKU4-1 | NC009019.1 
Bat SARS coronavirus HKU3-1 | DQ022305.2 
- Bat SARS-like coronavirus RsSHC014 | KC881005.1 
9 ^- Bat SARS-like coronavirus WIV1 | KF367457.1 
38«- SARS_SARS_coronavirus_|_NC004718.3 
- European turkey coronavirus 080385d | KR822424.1.1 

-Bulbul coronavirus HKU11-796 | FJ376620.1 


, Porcine deltacoronavirus isolate PDCoV/CHJXNI2/2015 | KR131621.1 


99 1— Porcine coronavirus HKU15 strain HKU15-44 I NC016990.2 
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Rhinolophus bat coronavirus HKU2 CH/GD-01/2017/P21 MF370205.1 
6C | • SADS-CoV/CN/GDTP/2017 S gene | MG605085 

• SADS-CoV/CN/GDGL/2017 S gene | MG775250 
Porcine enteric alphacoronavirus GDS04 | MF167434.1 

• SADS-CoV/CN/GDZW/2017 S gene | MG605086 

• SADS-CoV/CN/GDLS/2017 S gene | MG775252 

• SADS-CoV/CN/GDDL/2017 S gene | MG605084 
SADS-CoV FarmC | MF094683.1 
SADS-CoV FarmB | MF094682.1 
SADS-CoV FarmA | MF094681.1 
SADS-CoV FarmD | MF094684.1 

Bat coronavirus HKU2 | NC009988.1 
Bat coronavirus HKU2 strain HKU2/HK/46/2006 | EF203065.1 
Bat coronavirus HKU2 strain HKU2/HK/33/2006 | EF203067.1 
96 Bat coronavirus HKU2 strain HKU2/HK/298/2006 | EF203066.1 
BtRf-AlphaCoV/YN20121 KJ473808.1 


100 i BtRf-AlphaCoV/YN20121 NC028824.1 

- Bat coronavirus HKU9-1 | NC009021.1 

- Bovine coronavirus | NC003045.1 

- Human coronavirus OC43 strain SC2481 | KY983583.1 
r PHEV CoV USA-15TOSU1582 | KY419113.1 


£1 


100 l— Porcine hemagglutinating encephalomyelitis virus | NC007732.1 
Mouse hepatitis virus strain MHV-A59 Cl2 mutant | NC001846.1 
Murine hepatitis virus strain JHM complete genome | AC000192.1 
- Bat coronavirus HKU4-1 | NC009019.1 

-Bat coronavirus HKU5-11 NC009020.1 

-Coronavirus Neoromicia/PML-PHEI/RSA/20111 KC869678.4 

-Middle East respiratory syndrome coronavirus | NC019843.3 


-T 

100 1_r 

lOoL 


Bat SARS coronavirus HKU3-1 | DQ022305.2 
SARS coronavirus | NC004718.3 
Bat SARS-like coronavirus RsSHC014 | KC881005.1 


100 L Bat SARS-like coronavirus WIV1 | KF367457.1 
- European turkey coronavirus 080385d | KR822424.1.1 

100 r Porcinedeltacoronavirus isolate PDCoV/CHJXNI2/20151 KR131621.1 


L Porcine coronavirus HKU15 strain HKU15-44 | NC016990.2 

- Bulbul coronavirus HKU11-796 | FJ376620.1 

- Bat coronavirus HKU8 strain AFCD77 | EU420139.1 

- Rousettus bat coronavirus HKU10 | NC01887.1 

- Bat coronavirus 1A | NC010437.1 

Camel alphacoronavirus isolate Riyadh/Ry141/2015 | NC028752.1 
Human coronavirus 229E | NC002645.1 
229E-related bat coronavirus strain BtKY229E-1 | KY073747.1 


TL C 


- Human Coronavirus NL63 | NC005831.2 

— NL63-related bat coronavirus strain BtKYNL63-9a | NC032107.1 
100 I Transmissible gastroenteritis virus | AJ271965.2 

100 p JGEV Purdue P115 I DQ811788.1 
LpRCV ISU-1 | DQ811787.1 


- Bat coronavirus CDPHE15/USA/2006 | NC022103.1 
-Scotophilus bat coronavirus 512 | NC009657.1 


100 

97 


t 


Porcine epidemic diarrhea virus strain GDS01 | KM089829.1 
Porcine epidemic diarrhea virus strain CV777 | KT323979.1 
Porcine epidemic diarrhea virus| NC003436.1 


AplhaCoVI 


BetaCoV 


| GammaCoV 
DetalCoV 


AplhaCoV2 
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